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Abstract 

We investigate thermodynamics of the apparent horizon in f(R) gravity in the Palatini formalism 
with non-equilibrium and equilibrium descriptions. We demonstrate that it is more transparent to 
understand the horizon entropy in the equilibrium framework than that in the non-equilibrium one. 
£>> Furthermore, we show that the second law of thermodynamics can be explicitly verified in both 

phantom and non-phantom phases for the same temperature of the universe outside and inside the 
apparent horizon. 

"' PACS numbers: 04.50.Kd, 04.70.Dy, 95.36. +x, 98.80.-k 

u . 

> 

m 

(N 

*o 

in 

o 
o 



X 



* E-mail address: bamba@phys.nthu.edu.tw 
t E-mail address: gcng@phys.nthu.edu.tw 



I. INTRODUCTION 



There are two representative approaches to account for the current accelerated expan- 



sion 






2] of the universe |3Nl5|. One is the introduction of "dark energy" in the framework 
of general relativity. The other is the investigation of a modified gravitational theory, such 
as f{R) gravity 16Ml9|. In this study, we concentrate on thermodynamics in the modified 
gravitational theory The fundamental connection between gravitation and thermodynam- 
ics has been suggested by the discovery of black hole thermodynamics 20] with black hole 



entropy 2l| and Hawking temperature 22 



zon of de Sitter space has been explored 



ts application to the cosmological event hori- 
(for recent reviews, see 24|). It was shown 



that the Einstein equation can be derived from the proportionality of the entropy to the 
horizon area together with the Clausius relation in thermodynamics 



26r|3Q|. In Ref. |29[, i 



This consequence 



has been applied to various cosmological settings |26H30| . In Ref. [29|, it was shown that if 
the entropy of the apparent horizon in the Friedmann-Lemaitre- Robertson- Walker (FLRW) 
spacetime is proportional to the apparent horizon area, Friedmann equations follow from 
the first law of thermodynamics. 

Thermodynamics in dynamical black holes 
times 34J has also been developed. The first law of the ordinary equilibrium thermody- 



31 



331 ] and spherically symmetric space- 



namics has been discussed in f(R) gravity, scalar-tensor theories, the Gauss-Bonnet gravity 



and more general Lovelock gravity 35M39J]. while the second law has been studied in the 



accelerating universe, f(R) gravity, the Gauss-Bonnet gravity and the Lovelock gravity in 



Refs. 



40-42 



respectively. Thermodynamics in braneworld scenario [43|, loop quantum 
cosmology J44J, |45j and Hofava-Lifshitz gravity 46| as well as its properties of dark en- 
ergy 47H49I] and those in cosmological scenarios with interactions between dark energy and 



dark matter 



50 



51| have been a 



so examined. Incidentally, the horizon entropy in the four- 
dimensional modified gravity 52| and the quantum logarithmic correction to the expression 
of the horizon entropy in cosmological contexts 44|, I53m5| have been investigated. In the 
framework of general relativity, the considerations on the sign of entropy [56|, |57| and its 



evolution 



58| in the phantom phase have been performed. Moreover, the conditions 59] 



n 



that the black hole entropy can be positive in the f(R) gravity models [60h62| with passing 
the solar-system tests and cosmological bounds and the evolution of the black hole entropy 
in the ghost condensate scenario [63( have been studied. 



It was pointed out 64J that in f(R) gravity, a non-equilibrium thermodynamic treat- 
ment should be required in order to derive the corresponding gravitational field equation 



by using the procedure in Ref. 



25| . This point was reanalyzed in f(R) gravity 65| as well 



as scalar-tensor theories 66j and elaborated by considering the role of gravitational dis- 
sipation 67j. By taking into account the non-equilibrium thermodynamic treatment, the 



studies on the first [68 



and second 



691 ] laws of thermodynamics on the apparent horizon in 






generalized theories of gravitation and the investigations on thermodynamics |70j in a f(R) 
gravity model 7l| with realizing a crossing of the phantom divide from the non-phantom 
(quintessence) phase to the phantom one have been performed. The non-equilibrium cor- 



rection 



64J ] has been reinterpreted through the introduction of a mass-like function 72J and 



brmulation in Ref. 25| was also extended to the more 
, |77| (for related discussions, see 781] ). 



by other approaches 73M75I] . The 
general extended gravity theory [76 

Recently, it has been shown that it is possible to obtain a picture of equilibrium thermo- 
dynamics on the apparent horizon in the FLRW background for modified gravity theories 
with the Lagrangian density f(R,<j),X) (including f(R) gravity and scalar-tensor theories), 
where X = — (1/2) g^V ^(pV u (j) is the kinetic term of a scalar field <fi (V M is the covariant 
derivative operator associated with the metric tensor g^ u ), due to a suitable redefinition of 
an energy momentum tensor of the "dark" component that respects a local energy conser- 
vation 791 ] . 

The previous studies on thermodynamics in f(R) gravity have been executed in the 
standard metric formalism, in which the connection is defined in terms of the metric. There 
exist other formalisms to derive the gravitational field equation from the action, such as the 
Palatini formalism, in which the connection and the metric are assumed to be independent 
variables (for more detailed explanations on the Palatini formalism, see 10J). As a result, 
the gravitational field equation of f(R) gravity in t 



that in the metric one. Furthermore, f(R) gravity 80 



theories 



atini formalism is different from 



85 



84j and more general gravitational 



in the Palatini formalism have been examined. Bouncing cosmologies in 
the Palatini f(R) gravity [87J and the dynamical aspects of generalized Palatini theories of 
gravity 88[ have been also investigated. 

It should be mentioned that the Palatini f(R) gravity (with the exception of those with 
very heavily suppressed correction which would only become significant in the very far UV) is 



in conflict with Solar system tests 82j and particle physics 89j,|90j]. Moreover, it is plagued 



by surface singularities in stars and spherical matter configurations |91[ and it does not 
necessarily have a well posed Cauchy problem in the presence of matter 10] . However, it 
should be cautioned that some of the arguments may not be generic. For example, most of 
criticism were based on f(R) = R — 0/R n , where /3 and n are constants, while the against 
of the Palatini f(R) gravity based on particle physics may be premature as the Higgs sector 
is not well understood yet. 

On the other hand, although there have been many studies on thermodynamics in f(R) 
gravity in the literature, an equilibrium thermodynamics description of f(R) gravity in the 
Palatini formalism has not been done yet. In particular, it is interesting to ask whether 
the non-equilibrium and equilibrium descriptions of thermodynamics can be given in the 



Palatini formalism as the metric one 79j . For the interest of the formalism, in this paper 
we explore both non-equilibrium and equilibrium descriptions of thermodynamics in the 
Palatini formalism of f(R) gravity. We also show that the second law of thermodynamics 
in f(R) gravity can be explicitly verified in the phantom phase as well as the non-phantom 
(quintessence) one if the temperature of the universe inside the horizon is equal to that of 
the apparent horizon. We use units of &b = c = ft = 1 and denote the gravitational constant 
8ttG by k 2 = 8tt/M p1 2 with the Planck mass of M P1 = G~ l/2 = 1.2 x 10 19 GeV. 

The paper is organized as follows. In Sec. II, we explain f(R) gravity in the Palatini 
formalism and study the first and second laws of thermodynamics of the apparent horizon 
in f(R) gravity under a non-equilibrium description of thermodynamics. In Sec. Ill, we 
illustrate an equilibrium description of thermodynamics by redefining the energy density and 
pressure of dark components in f(R) gravity in the Palatini formalism. Finally, conclusions 
are given in Sec. IV. 

II. NON-EQUILIBRIUM DESCRIPTION OF THERMODYNAMICS IN f(R) 
GRAVITY IN THE PALATINI FORMALISM 

A. f(R) gravity in the Palatini formalism 

The action of f(R) gravity in the Palatini formalism with matter is written as 

'W) 



a xy/—g 



r\ 2 ^matter 



(2.1) 



where g is the determinant of the metric tensor g^ and £ m atter is the matter Lagrangian. In 
the Palatini formalism, the connection and the metric tensor g^ are treated as independent 
variables. We denote the Ricci tensor constructed with this independent connection as 1Z^ V 
and the corresponding Ricci scalar as TZ = g^ u TZ^ u . In general, TZ^ is different from the 
Ricci tensor R^ u constructed with the Levi-Civita connection of the metric. Here, the Ricci 
scalar R is given by R = g^R^ and f(TZ) is an arbitrary function of TZ. 

Taking the variation of the action in Eq. ( 12. ip with respect to g^, one obtains [lO] 

F{U)n, u - l -f{U)g, u = k 2 T^ , (2-2) 

where F(1Z) = df(lZ)/dlZ. Here, T^ a is the contribution to the energy-momentum 
tensor from all perfect fluids of ordinary matter (radiation and non-relativistic matter) with 
pi and Pf being the energy density and pressure of all ordinary matters, respectively. 

Varying the action in Eq. ( 12. ip with respect to the connection and using Eq. ( 12. 2ft . one 



finds [10] 



FG, V = K 2 T^ attcr) - \g»» (FK -f) + V,V U F - g^UF - JL 



V,FV U F - l -g^ (VF) 2 



(2.3) 

where G^ v = R^ u — (1/2) g^R is the Einstein tensor, V M is the covariant derivative operator 
associated with the Levi-Civita connection of the metric tensor g^, □ = g^V^V^ is the 
covariant d'Alembertian for a scalar field, and (VF) = g^V M FV U F . 

We assume the four-dimensional Friedmann-Lemaitre- Robertson- Walker (FLRW) space- 
time with the metric, 

ds 2 = h aP dx a dx p + f 2 dQ 2 , (2.4) 

where r = a(t)r, x° = t and x 1 = r with the two dimensional metric h a p = 
diag(— 1, a 2 (t)/[l — Kr 2 ]). Here, a(t) is the scale factor, K is the cosmic curvature, and 
dVt 2 is the metric of two-dimensional sphere with unit radius. In the FLRW background 
( 12 .4p . from Eq. ( 12. 3 p we obtain the following gravitational field equations: 

3F { R2 + ^) = n2pi + \( Fn -n- 3Hp - It ' (2 - 5) 

-2f(h-£\ =K 2 ( Pi + P i ) + F-HF-~ J (2.6) 



where H = a/a is the Hubble parameter and the dot denotes the time derivative of d/dt. 
We note that the perfect fluid satisfies the continuity equation 



Pi + 3H (p f + Pi) = 
Equations (12 .5p and (12. 6p can be rewritten as 



H 2 



K 



a z 



K 

3F 



(Pd + Pf) 



where pd and Pd are the energy density and pressure of "dark" components, given by 

1 3 F 2 

-(FK-f)-3HF-- — 

2 v j j 4 F 

1 - ■ 3 F 2 
— (FK - f) + F + 2HF 

2 v ; A F 







1 


IM 


~ 


K* 






1 


Pa 


= 


K* 



(2.7) 

(2.8) 
(2.9) 

(2.10) 
(2.11) 



Here, a hat denotes quantities in the non-equilibrium description of thermodynamics. Note 
that pa and Pd are originated from the energy-momentum tensor TJ^ U , defined by 



f(d) 

fJ,V 



\ 9lM , {FTZ -f) + V„V„F - g^UF 



IF 



V,FV U F - \ 9ilv (VF) 2 



where the Einstein equation is given by 



/~i /^(matter) _i_ rp(d) \ 



'fll/ 



Relation between TZ^ and R^ v as well as 1Z and R are given by [10 1 

1 



K 



t w 



V ^W,FW V F - I ( 



•t^ 



:^OF 



n = R+^- 2 (VF) 2 -^DF. 



(2.12) 
(2.13) 

(2.14) 
(2.15) 



In the FLRW background (12 .4p . the scalar curvature in the metric formalism is written as 
R — 6 ( 2/7 2 + H + iiT/a 2 1 . Using this expression and Eq. (12. 15ft . we get 



K = 6 [2H 



H + «) + itJ£ + P + 3HF 



(2.16) 



From Eqs. (I2.10p . (12. lip and (I2.16p . we find that pd and Pd satisfy the following equation: 



-^[H 2 + ^) 



(2.17) 



In f(R) gravity, since F^O, the right-hand side (r.h.s.) of Eq. ( 12 . 1 T[) does not vanish, so 
that the standard continuity equation does not hold. 



6 



B. First law of thermodynamics 

We study the thermodynamic property of f(R) gravity. The dynamical apparent horizon 
is determined by the relation h a ^d a fd^f = 0. According to the recent type la Supernovae 
data, it is suggested that in the accelerating universe the enveloping surface should be the 



apparent horizon rather than the event one from the thermodynamic point of view 40| . In 
the FLRW spacetime, the radius f A of the apparent horizon is given by 

f A ={H 2 + -j . (2.18) 

Taking the time derivative of this relation, we obtain 

Combing Eqs. (Q and (l2~T0|) . we find 

dr A = r A H(p d + P d + Pi + P^dt. (2.20) 



4vrG 
In general relativity, the Bekenstein-Hawking horizon entropy is expressed as S 



Aj (4G), where A = 4nr A is the area of the apparent horizon |20j-|22[. The Bekenstein- 
Hawking entropy S = Aj (4G) is a global geometric quantity which is proportional to the 
horizon area A with a constant coefficient 1/ (4G). This quantity is not directly affected by 
the difference of gravitational theories, i.e., the difference of the derivative of the Lagrangian 
density / with respect to R, F. On the other hand, in the context of modified gravity 
theories including f(R) grav ity a horizon entropy S associated with a Noether charge has 



been proposed by Wald 92[. The Wald entropy S is a local quantity defined in terms of 
quantities on the bifurcate Killing horizon. More specifically, it depends on the variation of 
the Lagrangian density of gravitational theories with respect to the Riemann tensor. This 
is equivalent to S = Aj (4G e ff), where G e ff = G/F is the effective gravitational coupling in 



f(R) gravity [93] 



By using Eq. (I2.20p and the Wald entropy in f(R) gravity in the Palatini formalism 94 1 



FA 

*=«r (2 - 21) 



we get 

-dS = 47rr A H (p d + P d + p { + P t ) dt + ^dF . (2.22) 

V / 2Gr 



2irr A 



We note that the Wald entropy in f(R) gravity in the metric formalism has the same form 



as in Eq. (12T2TJ) [92 



95 



The associated temperature of the apparent horizon has the following Hawking temper- 
ature T determined through the surface gravity n sg : 



T 



'SS 



r»,o 



2n 
1 



--9 a 



-hh af> d. 



J' 1 



2V-h 

where h is the determinant of the metric h a p and K sg is given by 



(2.23) 
(2.24) 



<sg 



•i-l 

ta 



2Hf A 



r 4(2H» + H+£ 



2irG„ 

-ta Pt - 3P T 



3F 



(2.25) 



with p t = Pd + Pi and P t = Pd + -Pf, denoted as the total energy density and pressure of the 
universe, respectively. It follows from Eq. ( I2.25P that if the total equation of state (EoS) 
w t = Pt/ Pt satisfies w t < 1/3, one has K sg < 0, which is the case for the standard cosmology. 
Using Eqs. ( 12. 23ft and (12.251) . the horizon temperature is given by 



T 



1 



ta 



2-kta \ 2Hf A . 
Multiplying the term 1 — f^/(2iJf^) for Eq. (I2.22|) . we obtain 



(2.26) 



TdS =-- 4irr A H ( p d + P d + p f + P { ) dt - 2nr A ( p d + P d + Pf + Pf ) df A + ^r\dF . (2.27) 



T 



In general relativity, the Misner-Sharp energy 
f(R) gravity, this may be extended to the form 



Refs. 



97 

68 



98 
69 



is defined as E = r A / (2G). In 



99 



1Q0|) 



72] (for related works, see also 



E 



? A F 
2G 



(2.28) 



From Eqs. ( 12TTH]) and (G^gjl . we find 

E 



f A F T/ 3P(P 2 + K/a 2 ) 

V- — — = V {pd + Pf) , 



(2.29) 



2G 8ttG 

where V = A7ir A /3 is the volume inside the apparent horizon. The last equality in Eq. ( I2.29P 
means that E corresponds to the total intrinsic energy. It is clear from Eq. (I2.29f) that E > 



8 



and therefore F > 0, which is consistent with the fact that the effective gravitational coupling 
in f(R) gravity G c s. = G/F should be positive 101 ] . In other words, the graviton is not a 
ghost in the sense of quantum theory [6l|. Using Eqs. ( 12. 7p and (j2.17p . we get 



dE = -4irr A H [p d + P d + p f + PA dt + 4<kt a (p d + p f ) dr A + —^dF 



I'A 



2G 



It follows from Eqs. (I2T27]) and (123D|1 that 



TdS = dE + 27ir A [p d + Pl - P d - P f )dr A + -^ (1 + ^f A T) dF . 



I'A 



2G 



By introducing the work density 

W 



32. 


33. 


66 



-- (r {mattcr)a Ph a p + f {d)a(3 h af} 



Eq. (I2.3ip is reduced to 



2 ( Pd + Pf - Pd - Pf 



I'A 



TdS = -dE + WdV + -£ (1 + 27rf A T) dF , 

2Gr 



which can be rewritten in the form 



(2.30) 



(2.31) 



(2.32) 
(2.33) 



(2.34) 



TdS + TdiS = -dE + WdV . 



(2.35) 



where 



di S = -lt± (1 + 2vrr^T) dF = - f | + 5 | ^ 



7T 



4F" 2 + F" + 3F/a 2 



G 



(F 2 + K/a?) ( 2F 2 + F + F/a^ 



t/F. 



(2.36) 



Equation ( I2.36P agrees with the result of Ref. [69( for K = obtained in f(R) gravity. 

The new term d^S can be interpreted as a term of entropy produced in the non-equilibrium 
thermodynamics. General relativity with F = constant leads to diS = 0, which implies that 
the first-law of equilibrium thermodynamics holds. On the other hand, in f(R) gravity the 
additional term in Eq. (12.36D appears because dF ^ 0. 



C. Second law of thermodynamics 

To investigate the second law of thermodynamics in f(R) gravity, we start with the Gibbs 
equation in terms of all matter and energy fluid, given by 

TdS t = d (p t V) + P t dV = Vdp t + (p t + P t ) dV , (2.37) 

where T and S t denote the temperature and entropy of total energy inside the horizon, 
respectively. Here, we have assumed the same temperature between the outside and inside 
of the apparent horizon. To obey the second law of thermodynamics in f(R) gravity, we 



require that [69J 



which leads to the condition [70J 

J= (h-—\ F>0. (2.39) 

In the FLRW background <®M . the effective EoS w cS is given by |8| w cS = -1-2H/ (3H 2 ). 
For H < 0, ty e ff > — 1, which corresponds to the non-phantom (quintessence) phase, while for 
H > 0, w c g < —1, which corresponds to the phantom phase. It follows from Eq. (I2.39J) that 
J > even in the phantom phase. As a consequence, the second law of thermodynamics in 
f{R) gravity can be satisfied in both phantom and non-phantom phases. This conseq uenc e 



is compatible with a phantom model with ordinary thermodynamics proposed in Ref. 101 ] 



III. EQUILIBRIUM DESCRIPTION OF THERMODYNAMICS IN f(R) GRAV- 
ITY IN THE PALATINI FORMALISM 

The reason why there exists a non-equilibrium entropy production term diS is that /)d 
and Pd defined in Eqs. f 1 2 . 1 j) and (12. lip obey Eq. (12.171) . whose r.h.s. does not vanish in 
f{R) gravity because F ^0. In other words, the standard continuity equation does not 
hold. In this section, we redefine the energy density and pressure of dark components to 
satisfy the continuity equation so that there is no extra entropy production term, which 
is referred as the equilibrium description. We will discuss such equilibrium description of 
thermodynamics in f(R) gravity in the Palatini formalism. 

10 



A. First law of thermodynamics in equilibrium description 



We rewrite Eqs. (J2.5P and ( 12. 6 j) in the following forms: 

FT 



3F F 



-2F ( H 



« (Pd + Pf) 



« 2 (p d + p d + Pf + p f ; 



(3.1) 
(3.2) 



where F is some constant, and p d and P d are the energy density and pressure of dark 
components redefined as 



Pd = — 



K Z 



(FK - /) - 3FF + 3 (F - F) [H 



K 



3_P 



■\ (FTZ-f) + F + 2HF - (F - F) f 



K 



2H + 3H z + —\ 



3F 2 



(3.3) 



(3.4) 



2 y ' \ a 2 J 1 F 

It follows from Eqs. ( 13 .3p and (13. 4p that the standard continuity equation can be satisfied 



p d + 3F (p d + F d ) = . 



In this representation, Eq. ( I2.20p becomes 

F 



4vrG 



df A = f\H (p d + F d + p f + P t ) dt . 



By introducing the horizon entropy S in the form 

F A 



S 



4G ' 



and using Eq. ( I3.6p . we obtain 



2irr A 



—dS = 4irr A H (p d + P d + pi + Pi) dt . 



From the horizon temperature in Eq. (I2.26P and Eq. (13. 8p . we find 

TdS = 4nr A H (p d + F d + p f + F f ) dt - 2nr 2 A (p d + P d + Pi + P { ) dr A 



By defining the Misner-Sharp energy as 

F r A 



E 



2G 



V (p d + Pf) 



we get 



dE = -4irf A H (p d + P d + p f + F f ) dt + 47rf£ (p d + p f ) df A 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



11 



where there does not exist any additional term proportional to dF on the r.h.s. due to 
the continuity equation ( 13. 5p . By combining Eqs. ( 13. 9 \ and (13. lip , we obtain the following 
equation corresponding to the first law of equilibrium thermodynamics: 

TdS=-dE + WdV, (3.12) 

where the work density W is defined by 

W = -{p d + pt-P d -P t ) . (3.13) 

As a result, an equilibrium description of thermodynamics can be derived by redefining the 
energy density p d and the pressure P d so that the continuity equation ( 13. 5p can be met. 

Note that the constant Fq can be chosen arbitrarily as long as Fq > in order to ensure 
that the sign of the Friedmann equation (13.11) does not change. It is considered that the 
natural choice is F = 1 because in this case the entropy S and the Misner-Sharp energy 
E reduce to the standard forms in the Einstein gravity: S = A/ (4G) and E = ?a/ (2G), 
respectively. 

From Eqs. (EH), (E2D and lEEBjl . we obtain 

~ F H (H -K/a 2 ) 

s = ^m\ (Pd + Pi + Pd + p f) = -V ^TW • (3 - 14) 

Thus, the horizon entropy increases as long as the null energy condition p t + P% = Pd + Pf + 
Pd + Pi > is satisfied. 

The reasons why the above equilibrium picture of thermodynamics can be realized are as 
follows: The first is that there exists an energy momentum tensor T^J satisfying the local 
conservation law V M T^ = 0. The second is that the entropy S is given by Eq. ( 13 ,7p as in 
general relativity, corresponding to the Einstein equation in the form 

^ = ^wr tcr) +^), (3.i5) 

where 

rW = ^{-\g,» [F (K-R) + F R -f} + V,V U F - g, u DF + (F - F) R, u 

(3.16) 



3 
2F 



V,FV V F - l 9ia> (VF) 2 



12 



It is clear that the local conservation of T^J follows from Eq. (J3.15P due to the relations 
V M G Ml / = and V^T^™ a er = 0. By using the conserved energy momentum tensor 7) 



/ ' (/ 



and the horizon entropy S defined in Eq. (13.71) and following the method in Ref. 25], the 
Einstein equation (13. 15ft can also be derived. 

It can be shown that the horizon entropy S in the equilibrium description has the following 



relation with S in the non-equilibrium description [791 ] : 

r^ dF _ 2n(F -F) H {i I -K/a 
2GT G (#2 + K/a 2 

By using the relations ( 12.361) and ( 13. 81) . Eq. ( 13.171) is rewritten to the following form: 



dS = dS + d, t S + T^dF "—^ , r „ , „, ox2 dt. (3.17) 



F F 4H 2 + H + 3K/a 2 K ' 

where 

* 6n AH 2 + H + 3K/a 2 dF 

^ = ~G H 2 + K,a 2 -R- (3 ' 19) 

The difference between S and S appears in f(R) gravity due to dF ^ 0, although S is 
identical to S in general relativity (F = F = 1). From Eq. (13. 18j) . we see that the change 
of the horizon entropy S in the equilibrium framework involves the information of both dS 
and diS in the non-equilibrium framework. 

In the flat FLRW spacetime (K = 0), the Bekenstein-Hawking entropy is simply propor- 
tional to the inverse squared of the expansion rate of the universe (S oc H~ 2 ) independent 
of gravitational theories. Hence S grows as long as H decreases, whereas the increase of 
H leads to the decrease of S (as in the case of superinflation). In other words, this prop- 
erty mimics the standard general relativistic picture with energy density pd and pressure P& 
of dark components defined in Eqs. ( 13. 3ft and ( 13.41) . It can be clearly understood that the 
Bekenstein-Hawking entropy grows for w t = Pt/pt > — 1, where p t = pd+pf and P t = P^+Pi, 
and that it decreases for w t < —1. 

The Wald entropy carries the information of gravitational theories through the depen- 
dence S oc FH~ 2 . For example, in a model of f(R) gravity f{R) = R + aR n , where a 
and n are constants, it follows that S oc H 2( - n ~ 2 ^ and hence S grows apart from n = 2 
because H increases (decreases) for n > 2 (n < 2). This evolution is different from that of 
the Bekenstein-Hawking entropy. The introduction of the entropy production term d^S in 
the non-equilibrium framework allows us to have a connection with the equilibrium picture 

13 



based on the Bekenstein-Hawking entropy as given in Eq. (j3.17p . Thus, it is considered 
that our equilibrium description of thermodynamics is useful not only to provide the general 
relativistic analogue of the horizon entropy irrespective of gravitational theories but also to 
understand the non-equilibrium thermodynamics deeper in connection with the standard 
equilibrium framework. 

Finally, we remark that the EoS of dark components in the non-equilibrium description 
is different from that in the equilibrium description. It follows from Eqs. (j2.10p and ( 12. lip 
that the EoS of dark components in the non-equilibrium description Wd = Pd/pd is given by 

2 (FK - /) - 12HF - 3F 2 /F 
On the other hand, by using Eqs. (13.31) and (j3.4p . the EoS of dark components in the 
equilibrium description w d = Pd/pd is expressed as 

-2 (FK - f)+AF + 8HF - 4 (F - F) ( 2H + 3H 2 + K/a 2 ) - 3F 2 /F 

w d = = ^ J- . (3.21) 

2 (FK - f) - Y1HF + 12 (F - F) (H 2 + K/a 2 ) - 3F 2 /F 

It is clear from Eqs. (I3.20p and (13.2ip that in general Wd ^ Wd except general relativity with 

F = F Q = constant. Thus, in order to compare the EoS of dark components predicted by 

a theory of f(R) gravity with the observations, it is necessary to derive the expression of 

the EoS of dark components in both non-equilibrium description and equilibrium one. This 

is another important physical motivation to consider both non-equilibrium and equilibrium 

descriptions. 

B. Second law of thermodynamics in equilibrium description 

To examine the second law of thermodynamics in the equilibrium description, we write 
the Gibbs equation in terms of all matter and energy fluid as 

TdS t = d ( Pt V) + P t dV = Vdp t + ( A + P t ) dV . (3.22) 

The second law of thermodynamics can be described by 

do dot ,„ , 

* + ^°- < 3 - 23 > 



which gives 



12nF H(H-K/a*) l 

> (3.24) 



G (H 2 + K/a 2 Y R 
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by using V = 4vr^/3, R = 6 (h + 2H 2 + K/a 2 \ and Eqs. ( I2~26j) . ([372]) and (l3~T4j) . It is 
clear that in the flat FLRW spacetime (K = 0), the second law of thermodynamics can be 
met in both non-phantom and phantom phases, which is the same as the non-equilibrium 
description. This also agrees with the argument in Refs. 47], |50(. We remark that the 
above consequence can be shown explicitly only for the same temperature of the universe 
outside and inside the apparent horizon. As a result, a unified understanding between 
non-equilibrium and equilibrium pictures of thermodynamics has been obtained. 

IV. CONCLUSION 

In the present paper, we have studied the first and second laws of thermodynamics of 
the apparent horizon in f(R) gravity in the Palatini formalism. We have explored both 
non-equilibrium and equilibrium descriptions of thermodynamics in f(R) gravity. The equi- 
librium framework is more transparent than the non-equilibrium one because in the equilib- 
rium description the horizon entropy is described by the single expression S, whereas in the 
non-equilibrium description it consists of two contributions, the change of the horizon en- 
tropy dS and the term produced in the non-equilibrium thermodynamics di§. We have also 
shown that the second law of thermodynamics can be satisfied in not only the non-phantom 
phase but also the phantom one, provided that the temperature of the universe inside the 
horizon is equal to that of the apparent horizon. 
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